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A  SAGE  LONG-RANGE-RADAB  SIMULATION  MODEL 

hy 

R.  L.  Dugas 

Ihere  have  been  several  treatises  on  pulsed  search  radar  (e.g.,  ^,7)^  and 
also  on  SAGE-type  radars  (e.g.,  5,  5)  >dilch  utilize  radar-data-processing 
machines.  'Ihe  problem  of  characterizing  the  output  data  is  virtually  impossi¬ 
ble  due  to  the  conqplex  factors  \dilch  Influence  the  propagation  of  radar  pulses 
and  their  reflection  from  aircraft  targets.  Often  (e.g.,  2,  3,  U,  6,  7)  a 
slnqpllfying  assusqptlon  is  made  vhlch  allovs  a  mathematical  treatment  of  the 
problem.  In  (6)  Swerling  treats  two  types  of  target  fluctuations:  non- fluctuat¬ 
ing  targets,  and  targets  with  Independent  Rayleigh-distributed  pulse-to-pulse 
fluctuations.  He  sxiggest  in  (6)  that  radar  data  for  other  types  of  fluctuations 
should  be  studied.  This  suggestion  is  one  of  the  objectives  in  developing  the 
Long- Range-Radar  (LRR)  simulation  model.  Because  of  the  difficulty  in  treating 
analytically  radar  targets  with  relatively  slow  scintillation  rates,  and  because 
of  the  complexity  of  the  SAGE  detection  logic,  it  was  decided  to  develop  a 
mathematical  model  using  random- saji^llng  techniques.  This  model  would  be  able 
to  generate  radar  data  for  virtually  any  type  of  target  scintillation.  The  out¬ 
put  of  the  model  would  be  sufficiently  detailed  so  that  all  properties  of  the 
data  could  be  studied. 

A  computer  program  designed  to  accang>li5h  these  objectives  has  been 

written  for  the  AN/fSQ-7  computer  at  Santa  Monica.  Ihe  model  simulates  data 

from  a  single  target  and  radar  as  the  target  follows  a  specified  flight  path. 

The  main  output  of  the  model  consists  of  sequences  of  quantized  video,  as 

they  would  be  generated  by  the  AN/fST-2  radar-data-processing  machine  when 

’^Throughout  this  document,  numbers  in  parentheses  indicate  items  in  the  list  of 
references  on  page  lU. 
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the  target  1b  scanned.  Quantized  video  Is  In  the  one  state  If  a  received 
radar  pulse  exceeds  a  threshold  In  the  AN/fST-2.  The  seimenees  of  qjuantlzed 
video  are  further  processed  "tiy  the  program  to  determine  If  the  target-detection 
threshold  has  been  attained.  If  so,  the  ccnputatlon  of  the  reported  target 
azimuth  Is  carried  out,  and  coo^pared  with  the  theoretically  correct  target 
azimuth  In  order  to  determine  the  azimuth- reporting  error.  The  statistics 
of  this  error  as  a  function  of  range  or  blip-scan  ratio  are  easily  determined. 
Die  reported  target  azimuth  Is  detezmlned  by  the  beam- splitting  technique 
used  In  the  AN/FST-2.  This  Involves  determining  the  azimuths  of  the  target's 
leading  and  trailing  edges  and  dividing  the  difference  iv  two.  The  difference 
Is  referred  to  as  run  length,  a  quantity  which  Is  transmitted  with  range  and 
azimuth  to  the  Direction  Center. 

The  model  congnites  run  length  for  each  detected  sequence  of  quantized 
video.  This  data  Is  useful  In  checking  the  properties  of  the  model  ou^ut 
since  there  Is  canq)arable  eniplrlcal  data  on  run- length  variations. 

The  flight  path  Is  divided  Into  range  Intervals  of  l6  n.m.  each.  This 
provides  sufficient  range  resolution  for  blip  scan  and  other  functions  of 
range.  With  a  radial  flight  path  and  a  target  speed  of  ^  knots,  about  one 
hour's  cooqniter  time  will  provide  enough  data  to  give  reasonably  high 
statistical  confidence  In  the  results.  This  time  estimate  assumes  that  both 
noimal  and  MTI  video  data  are  generated  for  all  ranges. 

The  model  simulates  cross-section  fluctuations  using  the  concept  of 

phase- sensitive  cross  section,  c^,  as  defined  In  (1).  This  Involves  q^eclfy- 

Ing  Individual  scatterers  on  the  target  and  their  movement.  If  o.  Is  the 
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cross  section  of  the  Jth  scatterer  on  the  target;  then 
a  =  i  y  vr  e^^j 

J 

vbere  the  Jth  phase  angle,  depends  on  the  position  of  the  Jth  scatterer, 
and  the  radar-vave  length,  X;  that  is 


\diere  d.  is  the  distance  from  the  Jth  scatterer  to  a  reference  plane  pexpendl- 
J 

cular  to  the  line  of  sight.  Of  course,  it  is  laipossible  to  specify  accurately 
the  cross  section  of  all  the  scatterers  for  all  aspects.  Ihe  target  cross 
section,  however,  can  be  modeled  at  L-band  wave  lengths  with  reasonable 
accuracy  by  considering  only  the  predominate  scatterers.  These  are  the  scatter¬ 
ers  whose  cross  sections  are  significant  fractions  of  the  average  target  cross 
section.  The  modeling  work  is  made  easier  by  assuming  that  two  or  more 
scatterers  in  close  proximity  can  be  represented  by  one  scatterer.  This  does 
not  significantly  change  the  fluctuation  rates  generated  by  the  model.  These 
criteria,  with  the  aid  of  other  techniques  of  (1),  cure  used  to  specify  the 
target  cross-section  input  data.  The  model  assumes  that,  although  the  cross 
section  of  a  scatterer  may  change  greatly  over  a  wide  aspect  angle  Interval, 
it  does  not  change  during  the  brief  e]g;>osure  time  \dien  the  radar  beam  sweeps 
past  the  target.  The  cross-section  fliustuatlons  during  e3g>osure  time  are 
generated  by  the  variations  in  the  dj  distances  as  a  result  of  the  taxget's 
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angular  movement  relative  to  the  reference  plane.  The  sharp  rise  in  cross 
section  at  certain  aspect  angles  due  to  the  edges  of  vlng  and  tall  surfaces  is 
accon5)lished  by  assuming  that  the  target  consists  of  only  two  nearly  equal 
scatterers  located  at  the  ends  of  each  edge. 

Hie  movement  of  the  scatterers  is  generated  entirely  by  yaw  movements. 
This  assur^tion  is  made  since  there  is  evidence  that  pitching^  rolling^  and 
air-frame  flexing  generally  contribute  little  to  cross-section  fluctuations 
centred  to  yaw  movements.  Ibe  yaw  movements  are  determined  by  two  Jaw  rates: 
the  yaw  rate  caused  by  the  movement  of  the  target  along  the  flight  path  and  a 
statistically  vaiying  yaw  rate  ^diich  takes  into  account  variations  due  to  air 
turbulence.  The  second  yaw  rate  varies  from  scan  to  scan  following  a  normal 
distribution.  Both  yaw  rates  remain  constant  during  exposure  time. 

Target  cross  section,  target  range,  and  radar  characteristics  are  com¬ 
bined  in  the  radar-range  equation  to  determine  predetector  signal-to-noise 
power  ratio.  The  following  decibel  form  of  the  radar-range  equation  is  used: 

SP/mP  =  a  +  2G  -  4r  +  K 

idiere  o  is  target  cross  section,  G  is  antenna  gain,  and  R  is  target  range,  all 
in  db  above  convenient  xmits.  The  antenna  gain  is  computed  from  the  expression 
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^ere  Is  the  antenna  gain  along  the  beam  axis,  W  Is  the  angle  between  the 
beam  axis  and  the  line  from  radar  to  target,  and  B  Is  a  constant  depending 
on  beam  width.  Ihe  values  of  G  for  values  of  W  corresponding  to  pulses  are 
precoo^ted  and  become  an  li^mt  table  for  the  model.  The  constant  K  Mcounts 
for  all  other  required  radar  characteristics  and  also  takes  Into  account  a 
imlt  conversion  factor.  K  contains  severed,  constants  idilch  are  difficult  to 
specllV;  however,  by  applying  the  model  to  some  situation  for  \dilch  a  blip* 
scan  curve  has  been  accurately  established.  It  Is  easy  to  euscurately  deter¬ 
mine  K. 

Predetector  slgnal-to-noise  ratio  (a)  Is  obtained  from 

■f »  2 

’>  SP/NP  =  10  log^  I 


Probability  of  pulse  detection  as  obtained  from  (?)  Is 
P(a)  vdv  e  — -g-  Ip(av) 


where  v  Is  the  nonnallzed  amplitude  threshold  in  the  quantizer  (the  pulse 
detector)  and  1^  Is  a  modified  Bessel  function  of  the  first  kind  and  zero 
order.  A  series  e:q>anslon  Is  used  to  calculate  P(a).  It  Is  given  by 


P(a)  =  1 


n=o  k=o 
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;Aiere  Pjj  is  the  probability  that  receiver  noise  alone  will  cause  the  pulse  de¬ 
tector  to  go  to  the  one  state.  Ihirteen  terms  of  the  first  sum  Insure  adequate 
accuracy  in  the  calculation  of  P(a).  In  SAGE  radars,  is  fixed  at  sojne  level. 
The  expression  for  P(a)  is  precon5)uted  and  tabulated  as  an  input  for  the  model. 
Singel-delay  moving-target  indicator  (MTl)  video  is  simulated  in  the  model  by 
multiplying  the  signal-to-noise  ratio  by  an  MTI-response  factor  before  deter¬ 
mining  the  probability  of  pulse  detection.  Die  response  factor  is  a  sinusoidal 
function  of  radial  velocity  (Vj^)  radar  prf  (P)  and  wave  length  (X).  It  is 
zero  at  blind  speeds  and  rises  to  a  maximum  of  O.9OO5  between  blind  speeds. 

The  expression  used  is 


The  model  determines  P(a)  for  every  pulse  of  every  scan  during  exposure  time. 
Before  and  after  exposure  time  it  is  assumed  that  P(a)  =  P^  for  enough  radar 
pulses  to  Insure  proper  detector  action.  Both  probabilities  are  sampled  to 
generate  sequences  of  quantized  video.  The  statistical  variation  in  the 
sequences  is  influenced  by  target-aspect  angle,  range,  scatter  cross  sections, 
scatter  movement,  euid  the  horizontal  antenna-gain  pattern. 

In  its  present  form,  the  model  is  only  two  dimensional.  It  appears  that 
introducing  a  third  dimension  in  describing  scatterer  movement  is  imnecessazy. 
The  height  effects  on  the  radar-gain  pattern  have  not  been  Included,  but 
could  be  taken  into  account  in  future  work.  The  most  serious  omission  is 
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the  multipath  effect;  hovever>  this  generally  occurs  only  vhen  the  terrain 
or  water  surrounding  the  radar  Is  smooth.  It  seems  reasonable  to  Interpret 
the  data  as  being  r^resentatlve  of  a  typical  situation.  Ihere  Is  evidence 
frem  the  results  of  an  application  that  bears  this  out. 

Ihe  model  has  been  exercised  for  one  type  of  radar  and  target.  Two 
flight  paths  have  been  used.  One  is  a  straight  radial  path  goixtg  from  maxi* 
mum  to  minimum  range.  The  other  path  consists  of  two  straight  paths  Joined  by 
a  circular  curve.  The  straight  paths  are  chosen  so  that  the  aspect-angle 
variation  covers  a  large  interval. 

Ihe  pulse-to-pulse  cross-section  fluctuations  were  quite  different  from 
those  of  either  (4)  or  (7).  Even  for  a  relatively  hl^  yaw- rate  standard 
deviation^  the  magnitude  of  cmsecutlve  cross  sections  could  not  be  considered 
ind^endent . 

Distribution  curves  for  crcss-sectlon  data  sanpled  from  the  pulse  cor¬ 
responding  to  the  center  of  the  beam  were  constructed.  For  the  radial  path, 
the  data  resembled  a  Poisson  distribution.  Using  the  other  path,  a  thirty- 
degree  change  in  aspect  angle  gave  data  which  was  more  like  a  Rayleigh 
distribution.  The  occurrence  of  Rayleigh-distributed  cross-section  data  la 
discussed  In  (9)  and  elsewhere.  In  (8)  It  is  pointed  out  that  L-band  cross 
sections  need  not  be  Rayleigh  distributed.  On  the  bMls  of  the  model  assump¬ 
tions,  and  because  of  the  lack  of  esplrlcal  L-band  cross-section  fluctuation 
data,  it  Is  reasonable  to  assume  that  the  data  generated  Is  realistic. 
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Blip- scan  curves  were  coo^lled  for  both  flight  paths.  data  for  the 
radial  path,  with  a  very  small  variation  in  aspect  angle,  agreed  well  with 
enplrical  blip- scan  curves.  ISie  path  with  a  large  variation  in  aspect  ai;gle 
produced  a  blip-scan  curve  which  diverged  considerably  frcm  that  of  the 
average  target.  This  data  is  shown  in  Figure  1.  Although  most  en^lrlcal 
blip-scan  curves  seem  to  be  compiled  from  radieLL- flight-path  data,  it  is  not 
unreasonable  to  find  that  the  blip- scan  function  changes  substantially  when 
the  flight  path  gives  a  wide  variation  in  aspect  emgle.  Generation  of  this 
kind  of  data  should  be  more  realistic  than  sampling  fixed  blip- scan  curves 
in  order  to  simulate  SAGE  data  trails. 

Ihe  distribution  of  run  lengths  for  the  radial- flight  path  agreed  well 
with  empirical  run- length  distributions  for  a  large  number  of  random  targets 
with  random  aspect  angles.  For  this  reason  it  was  decided  to  examine  this 
data  in  more  detail.  Ihe  conpoeltlon  of  the  sequences  was  further  analyzed 
by  coBiputing  distribution  curves  for  the  frequency  of  occurrence  of  various 
"solidities."  The  solidity  of  a  sequence  is  defined  as  the  maximum  number 
of  consecutive  hits  in  the  sequence.  Three  Ba]]g>le  spaces  were  used:  the 
space  consisting  of  all  data  for  all  ranges,  and  the  data  for  the  range  inter- 
veds  with  blip-scan  ratios  of  33  per  cent  and  53  per  cent.  Ihese  distribution 
curves  are  shown  in  Figure  2.  Although  the  beam  width  of  the  radar  corres¬ 
ponded  to  16  radar  pulses,  there  were  memy  sequences  with  solidities  greater 
than  this.  It  is  curiously  interesting  to  note  that  the  distribution  curves 
for  the  35  and  58  per  cent  blip-scan  range  intervals  sure  quite  simileo:.  Both 
extend  from  solidity  2  through  I?.  Both  have  median  values  which  are  between 


Figure  1 


Solidity  Distributions 


Figure  2 
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solidity  6  and  Tt  and  both  have  the  most  frequent  solidity  of  4.  Ihe  distri¬ 
bution  curve  for  the  space  of  all  the  data  Is  almost  uniform  between  solidities 
}  and  23.  !nils  shows  the  imlnltlatlve  nature  of  sequences  of  quantized  video, 
and  the  difficulty  In  characterizing  them.  To  say,  for  example,  that  a 
sequence  of  a  certain  number  of  consecutive  hits  would  be  representative  of 
a  certain  blip- scan  ratio  would  not  make  much  sense  on  the  basis  of  these 
results.  Another  procedure  is  to  say  that  the  hits  eure  uniformly  distributed 
over  the  width  of  the  beam.  Suppose,  for  example,  one  considered  the  space 
consisting  of  all  combinations  of  eight  hits  and  el£^t  misses,  nie  solidity- 
distribution  curve  for  this  space  Is  substantially  different  from  any  of 
those  of  Figure  2.  It  rises  to  a  sharp  peak  of  43.1  per  cent  at  solidity  3* 

The  radial-path  data  was  further  examined  to  determine  the  variation  In 
azimuth- reporting  Mcuinuiy.  It  was  found  that  the  mean  azimuth- rQ)ortlng 
error  varied  some\diat  with  received  signal  strength.  This  result  has  some 
Influence  on  the  best  azimuth-counter  preset.  The  best  preset  Is  discussed 
In  (2),  \diere  this  dependence  Is  not  noted.  An  axialysls  of  the  radlal-fli^t- 
path  data  showed  that  the  mean  azimuth- r^orting  error  became  about  one-half 
azimuth  count  Icurger  as  the  blip- scan  ratio  varied  from  Its  highest  to  lowest 
values.  This  means  that,  because  of  signal- strength  variations,  random 
azimuth-reporting  errors  can  Introduce  a  one-azimuth-count  bias,  even  when 
the  best  azimuth  counter  preset  is  used. 

The  standard  deviation  of  azimuth- r^ortlng  errors  vas  conputed  for  the 
radial- flight  path  as  a  function  of  range.  For  eonparatlve  purposes,  the 
model  was  rerun,  using  a  target  with  a  single  scatterer  with  the  same  cross 
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section  as  the  average  target  cross  section  of  the  previous  runs,  and  the 
standard  deviation  of  azlouth> reporting  errors  vas  congnited.  This  gave  data 
for  a  conparable  nonfluctuating  target.  These  results  are  shown  on 
Figure  3  idilch  also  includes  plots  of  the  theoretical  mlnlmun  standard  devi¬ 
ation  (6)  for  the  same  target  and  radar,  assuming  that  the  target  Is  either 
not  fluctuating  or  has  Independent,  Raylelgh-dlstrlbuted,  pulse-to-pulse 
fluctuations.  It  Is  Interesting  to  note  that  both  sets  of  data  are  similar, 
except  for  extreme  ranges.  The  difference  In  the  levels  of  the  two  sets  of 
data  could  probably  be  eiqplalned  by  quantlzatlco  errors  and  the  fact  that  the 
AN/FST-2  usee  a  detector  logic  ihlch  Is  not  quite  as  good  as  the  Ideal 
observer  theory  used  In  (6).  The  difference  In  the  shapes  of  the  curves  at 
extreme  ranges  Is  difficult  to  e]g>laln.  The  model  data  for  the  fluctuating 
target  does  not  continue  to  rise  at  extreme  ranges  as  does  the  theoretical 
data,  but  tends  to  converge  vith  the  nonfluctuatlng  data.  It  should  be 
pointed  out  that  the  sasq>le  size  gets  smeCLler  with  range;  but  an  additional 
set  of  runs  for  the  fluctuating  target  produced  the  same  effect.  More  study 
of  the  model's  behavior  Is  needed  to  e]Q>lain  this  effect. 

From  (2)  a  coaq>arable  estimate  of  the  standard  deviation  of  the  azimuth- 
reporting  error  for  a  "weak  tcurget"  has  been  obtained.  The  standard  deviation 
is  equivalent  to  0.0^  beam  widths.  This  value  seems  to  be  low  compared  to 
the  rest  of  the  data  of  Figure  To  the  author's  knowledge  there  Is  no 
empirical  data  which  could  be  compared  with  Figure  3* 

Scan-to-scan  data  for  the  radial  path  with  both  the  fluctuating  and 
nonfluctuating  target  was  obtained  by  the  model.  The  data  was  reduced  to 
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detemine  the  scan- to- scan  persistency  of  detectable  8eq.uenees  as  the  blip- 
scan  ratio  changes.  It  has  been  shown  (8)  that  sane  enplrlced  data  exhibits 
scan- to- scan  persistency  idilch  can  be  characterized  as  a  second-order  Markov 
process.  Ihe  expression  describing  the  process  Is: 

-  1  -  a  (1  -  p) 

where  >  probability  that  a  detectable  sequence  is  followed  by  a 
detectable  sequence. 

P  =  bllp-scan  ratio 
a  -  a  constant 

Different  values  of  a  produce  stral^t  lines  on  Figure  4^  above  the 
diagonal  and  passing  through  ^  1>  P  =  1.  Data  along  the  diagonal  (a  s  i) 
Indicates  scan- to- scan  ind^endence  (l.e.,  =  p)> 

Figure  4  shows  that  the  model  data  exhibits  no  scan-to-scan  correlation 
In  addition  to  that  liqplled  by  p.  None  was  to  be  e:q)ected  for  the  non- 
fluctuating  target.  Hie  fluctuating  target  produced  variations  idilch  were 
high  enough  to  be  considered  Ind^endent  fren  scan  to  scan.  It  Is  the 
opinion  of  the  author  that  a  small  yaw-rate  standard  deviation,  idien  used 
with  a  flight  path  having  a  very  low  rate  of  change  In  aspect  angle  may  give 
rise  to  scan-to-scan  data  which  lies  well  above  the  diagonal  of  Figure  4. 

Hils  and  other  aspects  of  the  model's  behavior  remain  to  be  studied  In 
further  applications  of  the  model. 
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Halar  Range  Ccegiutera 

Rqitorti  that  a  aathenatlccd  model 
using  randoa-eaapllng  techniques  has 
been  developed,  to  slaqilify  the 


analjrtleal  treatment  of  radar  tcurgets 
with  relatively  slow  scintillation 
rates.  Ihls  was  necessaxy  because  of 
the  coiiq>lezl'ty  of  the  SAQE  detection 
logic.  States  that  thlb  model  is 
able  to  generate  radar  data  for 
virtually  any  type  of  target 
scintlUatlon.  Describes  a  ecopu'ter 
program  written  for  the  AR/^SQ-7 
conqni'^^i^  at  Santa  Monica. 
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